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Any quantum system, such as those used in quantum information, magnetic resonance, or the like, 
is subject to random phase errors that can dramatically impact the fidelity of a desired quantum 
operation or measurement [T] . In the context of quantum information, quantum error correction 
techniques have been developed to correct these errors, but resource requirements are extraordinary. 
Realizing a physically tractable quantum information system will thus benefit significantly if qubit 
error rates are far below the so-called Fault-Tolerance error threshold pQ, predicted to be of order 
10~ 3 — 10 -6 . The need to realize such low error rates motivates a search for alternate strategies to 
suppress errors in quantum systems|i2]. We present experimental measurements on a model quantum 
system that demonstrate our ability to dramatically suppress qubit error rates by the application 
of optimized dynamical decoupling [3, 4 5, 6, 7 8, 9 pulse sequences in a variety of experimentally 
relevant noise environments. We provide the first demonstration of an analytically derived pulse 
sequence developed by Uhrig [10], and find novel sequences through active, real-time experimental 
feedback. These new sequences are specially tailored to maximize error suppression without the 
need for a priori knowledge of the ambient noise environment. We compare these sequences against 
the Uhrig sequence, and the well established CPMG-style spin echo |lll 112] , demonstrating that our 
locally optimized pulse sequences outperform all others under test. Numerical simulations show 
that our locally optimized pulse sequences are capable of suppressing errors by orders of magnitude 
over other existing sequences. Our work includes the extension of a treatment to predict qubit 
decoherence |13l I14| under realistic conditions, including the use of finite-duration, square ir pulses, 
yielding strong agreement between experimental data and theory for arbitrary pulse sequences. 
These results demonstrate the robustness of qubit memory error suppression through dynamical 
decoupling techniques across a variety of qubit technologies 12l [T5l 1161 1171 119] . 



We consider classical phase randomization of our qubit 
due to the action of the environment as the dominant 
source of errors, a situation germane to many quantum 
systems. Accordingly, we may write a Hamiltonian as 
H = + (3(t)]az, where ft is the qubit splitting, /3 
is a classical random variable (T4J [18], and <Jz is a Pauli 
operator. As in reference |14j . we may write the time 
evolution of a superposition state initially oriented along 
Y under the influence of this Hamiltonian as |^(t)) = 
^ (e -*m/2 e -i />(*')*' | T) + e int/2j f* WW ^ 

with | 1) and | J.) the qubit states. The term (3(t) adds 
a random phase to the time evolution of the state, cor- 
responding to a random rotation around the equator of 
the Bloch sphere. The buildup of such a random phase 
results in decoherence as an observer loses track of the 
Bloch vector in the equatorial plane. However, the appli- 
cation of a 7r pulse around X (henceforth denoted ttx), 
at time t' — t/2 will approximately time-reverse the state 
evolution, so long as fluctuations in [3 are slow relative to 
the allowed free-precession time of the qubit. This is the 
basis of the Hahn spin echo, a fundamental technique for 
preserving coherence in nuclear magnetic resonance and 
electron spin resonance systems 

Given an arbitrary noise power spectrum Sa{uS), we 
would expect that the Hahn echo acts as a high-pass fil- 
ter, mitigating phase errors associated with slowly vary- 
ing Fourier components of f3. It was shown that this 
general interpretation can be extended to multipulse se- 
quences. Following references [13] and [T3], for any n 



pulse sequence of total qubit evolution time r, the coher- 
ence of the state is given as W(t) — |(ov)(£)| = e~ xtyt \ 
where angled brackets indicate an expectation value, and 
the overline indicates an ensemble average in a frame ro- 
tating with frequency f2. In this expression, 

oo 
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where the filter function F(lut) contains all information 
about how the pulse sequence will preserve qubit coher- 
ence under the influence of Sp{uS). F(lot) is calculated 
from F{ujt) = \y n (ijjT)\ 2 , where y n (ujT) is the Fourier 
transform of the time-domain filter function, y n (t) (Figs, 
lb & lc). The time domain filter function alternates be- 
tween +1 and -1 for successive free precession periods. 

We begin by studying two distinct pulse sequences 
as will be described in order below: CPMG, and UDD 
(pulse spacings illustrated in Fig. la). CPMG is an 
extension of the Hahn spin echo to a multipulse form, in- 
corporating evenly spaced tt pulses about an axis rotated 
90 degrees from the direction imparting the initial (§)x) 
(Fig Id). This sequence has been shown to be robust 
against a variety of phase and rotation errors, and does a 
particularly good job at refocussing the Bloch vector [3J. 

Uhrig discovered [TU] that for an n pulse sequence, it 
is possible to modify the sharpness of the filter function, 
and hence the efficiency of noise suppression relative to 
CPMG, simply by changing the relative positions of the 
pulses within the sequence (Figs, la & lc). He developed 
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FIG. 1: CPMG and UDD pulse sequence schematics, (a) 
Fractional pulse locations, Sj, of CPMG (solid markers) and 
UDD (open markers) sequences as a function of pulse number, 
n. (b) Examples of the time-domain filter function, y§(t), for 
the CPMG and UDD pulse sequences with 6 7r-pulses. The 
dotted line represents the time-domain filter function assum- 
ing delta-function n pulses, while the solid line represents the 
time-domain filter function accounting for a finite 7V. (c) Log- 
log plot of the log of the filter function, F(ujt) = |£/6(o;t)| , 
for sequence length r =30 ms, and tv =185 /xs. Dotted lines 
indicate filter function with delta-function tt pulses, solid lines 
account for finite t„-. (d) Generalized schematic of an experi- 
mental sequence showing key procedures and designations of 
key times. 



an optimized sequence for the spin boson model (appro- 
priate for many solid-state systems), UDD, designed to 
increase the suppression of errors at short times — the 
so-called "high-fidelity" regime (the UDD sequence was 
later shown to be general [DTI H71 US]). Pulses are no 
longer evenly spaced, and deviation from the CPMG se- 
quence grows with n (CPMG and UDD are equivalent 
for n=2). 

The experimental noise spectrum can vary signifi- 
cantly between qubit implementations (e.g., semicon- 
ducting quantum dots vs. ultracold atoms), and as such, 
it is not sufficient to find and test a single pulse se- 
quence that suppresses noise for only one characteristic 
spectrum. For example, in superconducting qubit sys- 
tems [TJl HH 122] i fluctuating electric charges and 



spin centers produce noise spectra varying as 1/w. By 
contrast, a spin-boson model for quantum noise in semi- 
conducting quantum dots suggests the presence of noise 
with an ohmic spectrum, Sg(cj) cx u>, and a sharp 
cutoff [151 1221 125 • In order to test the efficacy of any 
pulse sequence, one must develop a method and testbed 
capable of exhibiting a variety of realistic noise environ- 
ments. In this manner we may think of the testbed as 
being a model quantum memory capable of simulating 
nearly any other. 

We realize such a model quantum memory in an array 
of ~ 1000 9 Be + ions in a Penning ion trap 25j. Previous 
experiments have demonstrated that under appropriate 
conditions these ions form 2D or 3D arrays with well de- 
fined crystal structure [Ml 123 EB] (Inset, Fig. 2a). This 
system mimics dense qubit arrays that might be con- 
structed in a variety of physical incarnations including 
atomic systems, superconducting circuits, and semicon- 
ducting quantum dots. 

The qubit states are realized using a ground-state 
electron-spin-flip transition as shown in Fig. 2a (also 
see Methods). Qubit operations are achieved by directly 
driving this ~124 GHz transition via a quasi-optical mi- 
crowave system (Fig. 2b). We drive cry rotations by 
the application of square microwave pulses to the qubit 
(Figs. Id, 2c, & 2d), while controlled oz rotations are 
achieved by detuning the microwaves from resonance and 
allowing the qubit state to freely precess for a precise 
time period (Fig. 2e). These methods are particularly 
well-suited to dynamical decoupling studies compared 
to laser-mediated qubit rotations due to the absence of 
spontaneous emission in the microwave/millimeter-wave 
regime. Doppler cooling of ion motion along the axis of 
the Penning trap [33] using UV laser light red-detuned 
from an atomic transition yields ion temperatures of or- 
der 1 mK. State initialization occurs via optical pumping 
(Fig. 2a), and state readout is achieved by fluorescence 
detection on the same cycling transition used for cooling 
[3"U] . We are able to initialize the system in a pure state 
with high fidelity, and perform a strong projective mea- 
surement, despite the use of ensemble techniques [3 1. 

We apply pulse sequences consisting of a few to more 
than 1000 ttx pulses. We have successfully extended our 
qubit coherence time (i.e. 1/e decay time) from approx- 
imately 1 ms as measured via Ramsey free-induction de- 
cay, to over 200 ms using 500 irx pulses in a CPMG 
sequence. We find that for large n, our minimum mea- 
sured deviation from full projection to | |) is of order 10 
% (likely due to long-time amplitude instability of our 
microwave system). Accordingly, we focus primarily on 
sequences with n < 10, allowing us to compare pulse se- 
quences in a regime where the minimum error rate is < 1 
%. 

We apply the CPMG and UDD sequences for various 
pulse numbers (Figs, la & Id), and measure ensemble- 
averaged state decoherence due to ambient magnetic field 
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FIG. 2: 9 Be + qubit structure and coherent control, (a) Rel- 
evant atomic structure of 9 Be + at 4.5 T. Qubit states labeled 
with (m_f,mj), the nuclear and total spin projections of the 
atom along the quantization axis. (Inset) Strobed optical im- 
age of ion fluorescence showing hexagonal-close-packed order 
with spacing ~10 (tin. (b) Schematic block diagram of the mi- 
crowave system used to drive qubit rotations. PLL = Phase 
Lock Loop, DRO = Dielectric Resonator Oscillator, CPL 
= Coupler. 77 MHz PLL reference provided by computer- 
controlled DDS or the beat note of DDS and a frequency 
modulated synthesizer during noise injection (see Methods), 
(c) Qubit transition driven via a square microwave ir pulse 
of Ttt « 185 (is. Readout consists of fluorescence detection 
immediately before and after application of microwave pulses 
followed by normalization to account for slow fluctuations in 
laser intensity (Normalized Counts). Each data point consists 
of the average of 20 experiments, (d) Rabi oscillations driven 
on-resonance give a decay time of 30-40 ms. tv is tunable via 
an input attenuator, (e) Ramsey fringes measured for two dif- 
ferent microwave detunings, 50 kHz and 100 kHz (tt/2 pulses 
were driven on-resonance and the microwaves detuned dur- 
ing the free- precession period). Detuning was controlled by 
switching between pre-programmed frequency output profiles 
via a TTL input on the DDS. 



fluctuations as a function of total free-precession time, 
as shown in Fig. 3a (see Methods). The ambient mag- 
netic field noise in our high-field superconducting mag- 
net is measured directly, varying approximately as 1 jui 2 
(Sp(uj) oc 1/w 4 ) with additional sharp spurs of undeter- 
mined origin, including a prominent feature at ~153 Hz. 
CPMG appears to outperform the analytically optimized 
UDD sequence for all n < 10 in this noise environment, 
which has a soft high-frequency cutoff [13] . Further, the 
data in Fig. 3a demonstrate that it is possible to extend 
the qubit coherence time by adding n pulses, as expected. 

Fitting our experimental data requires that we account 
for finite 7r-pulse durations in expressions for the pulse 



sequence filter function, diverging from the zero-pulse- 
length assumptions made in most literature on dynami- 
cal decoupling [HJ [19j [32] ■ We assume that dephasing is 
negligible during the application of a ttx pulse, based 
on the observation of Rabi-flopping decay times more 
than an order of magnitude longer than Ramsey free- 
induction decay times, and build on the theoretical de- 
scriptions of Uhrig[13 and Cvwinski[14j. The above as- 
sumption leads to the insertion of a delay, t„, between 
each free-precession time, during which the filter function 
in the time domain has value zero (rather than ±1 as in 
[TcH IT4l l3"2"]). Moving to the frequency domain, we may 
write the filter function of an arbitrary n-pulse sequence 
as 

F{WT) = \y n {uT)\ 2 

= |1 + (-l) n+1 e iuT (2) 

n 

+ 2 £(-1)V^ t cos(<^cjt/2)| 2 , 
i=i 

where SjT is the time of the center of the j th nx pulse, r 
is the sum of the total free-precession time and 7r-pulse 
times, and (j)^ is the ratio of tv compared to the total 
sequence length r (Figs, lb, lc & Id). To this order of 
approximation, all information pertaining to finite pulse 
lengths is accounted for by the simple addition of a co- 
sine term in the equation. 

Fits to experimental data show good agreement with 
theory. In Fig. 3a, the free fit parameters are the overall 
noise strength and the relative strength of the 153 Hz 
spur in our noise spectrum. This spur can be observed 
to slowly change amplitude in real-time, and is entirely 
responsible for the plateau-like feature we see in our de- 
coherence curves. Increasing the strength of this spectral 
feature changes the plateau-like feature to a rounded hill 
of increasing height. We believe that deviations between 
our experimental data and fitting functions are domi- 
nated by slow as well as discontinuous changes in the 
ambient noise environment. 

Numerical simulations suggest that in the "high- 
fidelity regime," UDD is capable of significantly outper- 
forming CPMG even in this noise environment, suppress- 
ing errors by several orders of magnitude to yield ultimate 
fidelities in excess of 99.99 % [HUE]. In order to empha- 
size the differences between these sequences, and over- 
come limitations in measurement fidelity, we artificially 
inject noise to simulate systems where UDD outperforms 
CPMG in the lower-fidelity regime (see Methods). 

We inject noise with ohmic and 1/cj power spectra 
and test the relative performance of UDD and CPMG 
(Figs. 3b & 3c). Our data indicate that the UDD se- 
quence dramatically outperforms CPMG in the presence 
of noise with an ohmic spectrum and a sharp cutoff [10 
- a significant departure from results under ambient 
noise. By contrast, using a 1/w power spectrum with 
a sharp cutoff, we find that over the entire range of ac- 
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FIG. 3: Pulse sequence performance in the presence of var- 
ious noise spectra, (a) Decoherence traces as a function of 
free-precession time for CPMG and UDD and ambient mag- 
netic field noise. Phase errors are manifested as nonzero flu- 
orescence detection. Traces saturate at 0.5 for total phase 
randomization. Traces for n=4,5 offset by 0.5 units for clar- 
ity. Each data point corresponds to 50 averages. Fits utilize 
ambient noise spectrum (see inset) and two free parameters - 
overall noise scaling (a) and scaling factor for the noise spur at 
~153 Hz (7). Extracted fit parameters yield a = 1.46 ± 0.31 
and 7 = 0.23 ± 0.04. Individual fit uncertainty is ~0.05- 
0.07 for each trace. Variations in extracted noise intensities 
are consistent with direct observation of slow fluctuations in 
the noise spectrum, including significant changes in the am- 
plitude of the ~153 Hz spur. (Inset) Ambient noise spec- 
trum measured via a solenoid embedded in our NMR magnet, 
(b) Performance of UDD and CPMG for Ohmic spectrum 
with 500 Hz cutoff (see inset). UDD outperforms CPMG 
by increasing factors as the noise intensity is increased. Sin- 
gle fit parameter: noise scaling, a. Extracted noise scaling: 
aohmic = 1-49 ± 0.19, expected a =1.65 for this noise inten- 
sity, (c) Sequence performance in a 1/lo noise spectrum with 
a sharp cutoff. Extracted noise scaling: aii/ u = 3.32 ± 0.15. 
Insets to (b) and (c)) Smoothed, measured noise on a log-log 
plot, measured using Vn = 0.7 V. Red line represents in- 
tended noise envelope up to high-frequency cutoff. Measured 
noise spectra extending from .01 Hz to 1 MHz used for fitting, 
but truncated for visibility. 



cessible noise intensities, CPMG performs similarly to 
UDD. As expected, peak noise power at low-frequencies 
is well-filtered by both sequences, yielding longer over- 
all coherence times for the 1/lo spectrum relative to the 
application of an Ohmic spectrum, and consistent with 
theoretical work [13]. Fits to theoretical expressions for 
the UDD and CPMG filter functions incorporating the 
effects of finite r w show strong agreement using a single 
fitting parameter: a, a scaling factor for the overall noise 
level. 

These experimental measurements are reproducible be- 
tween different ion clouds, and do not appear to vary sig- 
nificantly with changes in ion number. Additionally, the 
UDD sequence appears robust to rounding errors associ- 
ated with our pulse timing accuracy of 50 ns, and per- 
forms similarly for all input states on the equator of the 
Bloch sphere (not shown). We often used evenly spaced 
pulses (Periodic Dynamical Decoupling [4]) as a bench- 
mark; in all cases it performs worst of the sequences un- 
der consideration. Finally, in all noise environments and 
for all pulse sequences, neglecting the effects of finite T n 
dramatically diminishes our fit quality and leads to wide 
variations in the extracted noise strength, especially in 
regimes where the free-precession period is comparable 
to or less than t w . 

We extend the ideas of Uhrig by developing novel Lo- 
cally Optimized Dynamical Decoupling (LODD) pulse se- 
quences that are tailored to a given experimental noise 
environment and a fixed free precession time. We employ 
the Nelder-Mead simplex method for optimization in an 
n-dimensional space, manipulating the relative pulse po- 
sitions in the sequence. In our numerical simulations the 
Nelder-Mead method converges to an optimized pulse se- 
quence for given initial conditions in less than 100 iter- 
ations, providing a suite of locally optimized pulse se- 
quences (Inset, Fig. 4a). 

The efficacy of any optimization method in determin- 
ing LODD pulse sequences, however, is limited by the 
degree to which Sp(co) is known. In addition to gen- 
eral noise models, the specific measurement hardware 
and laboratory environment will add frequency compo- 
nents to Sp(uj) that are not easily measured or predicted 
by theory. We therefore modify the Nelder-Mead search 
algorithm by replacing numerical calculation inputs with 
the results of experimental measurements, resulting in a 
form of real-time, active, experimental feedback for quan- 
tum control [55] . This extension allows us to find opti- 
mized pulse sequences without definite knowledge of the 
experimental noise environment. 

Experimental optimization results for an Ohmic spec- 
trum are displayed in Fig. 4a, along with a schematic 
depiction of the resultant optimized pulse sequence com- 
pared with UDD (inset). We find that small modifica- 
tions in the pulse positions (for both ohmic and 1/lo spec- 
tra) produce significant improvements for the selected 
total free-precession time, and fits to theoretical expres- 
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FIG. 4: Nelder Mead pulse-sequence optimization using the 
ohmic spectrum in Fig. 3b. (a) Experimental measure- 
ments of 67r-pulse UDD and experimentally optimized pulse 
sequence derived from UDD on a logarithmic scale. Fit pa- 
rameters: ajjDD = 0.69, (XExp.Opt. = 0.77. Green arrow 
indicates the total sequence length for which experimental 
optimization was performed. The measured error rate us- 
ing the LODD sequence is ~ 7x lower than that obtained 
using CPMG. Blue dotted line indicates numerical simula- 
tion of decoherence with the same measured noise spectrum, 
a — 0.77, and with the pulse sequence optimized for each 
total sequence length, r. (Inset) Optimized sequence as a 
function of total free-precession time (solid markers). Black 
X's indicate the n pulse positions of the single experimentally 
optimized sequence. Ticks and horizontal grids correspond 
to original UDD pulse positions, (b) Filtration benefits of 
experimentally optimized pulse sequence tailored to injected 
noise environment. Noise power spectrum (upper trace) and 
logarithmic improvement of optimized filter function relative 
to UDD filter function (lower trace). Positive numbers (also 
background shading) indicate regions of the noise spectrum 
where the optimized sequence provides noise suppression su- 
perior to UDD. 



sions for qubit decoherence yield strong agreement. In 
the data presented we are able to suppress the qubit er- 
ror rate by nearly an order of magnitude over CPMG 
by the use of a LODD pulse sequence (a factor of five 
relative to the starting sequence, UDD). We also show 
via numerical simulation that it is possible to continu- 
ously vary the pulse sequences such that the filter func- 
tion is optimized for all total free precession times (blue 



dashed line, Fig. 4a), yielding sequences which outper- 
form CPMG and UDD in the high-fidelity regime by an 
order of magnitude. 

In summary, we have demonstrated the efficacy of 
phase error suppression via optimized dynamical decou- 
pling pulse sequences applied to a model quantum mem- 
ory. This experimental system has been employed to test 
various pulse sequences under experimentally realistic 
noise environments, yielding good agreement with theo- 
retical predictions for qubit decoherence under sequence 
application. Finally, we have developed a real-time ac- 
tive feedback technique to experimentally produce lo- 
cally optimized pulse sequences without any knowledge of 
the experimental noise environment, making these meth- 
ods generally useful to any experimental setting. Future 
studies will incorporate the concatenation of optimized 
pulse sequences [Sj [53] , and will examine the use of pulse- 
shaping for further error reduction. Our results provide 
key ingredients of a quantum toolkit which will make the 
production of a functional and useful quantum computer 
more realistic. 

METHODS 
Experimental control 

In our experiments we control several parameters 
using a centralized programmable- logic-device (PLD)- 
based pulse generator. We have pulsed control over de- 
tection/cooling and repump lasers, pulse counters, mi- 
crowave reference frequency, and a microwave shutter 
which provides ^25 dB of isolation in the off state. All 
relevant pulse sequences are generated in software on a 
controlling PC and sent to the PLD. Similarly, various 
settings are sent directly to the DDS which serves as a ref- 
erence for our microwave system, and to the data acqui- 
sition system which records fluorescence detection events 
from a phototube. An operator command launches the 
experiment which is synchronized to main line frequency 
and triggers the PLD to begin its pulse sequence. A 
complete experiment including state initialization, mea- 
surement, and the application of the pulse sequence it- 
self takes approximately 200ms. The experiment is con- 
ducted a fixed number of times and the results recorded 
and averaged before new experimental conditions are up- 
loaded. 

Pulse sequence application 

An experiment begins with initialization of the qubit 
in | t) (Fig. Id). The qubits are then rotated to the 
equator using a ^ pulse around the X axis, (5)x> form- 
ing a superposition state along Y. The microwaves are 
detuned during the interpulse period, and the qubits are 
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allowed to precess about Z at a rate controlled by the 
detuning from resonance (see Figs. Id & 2e). Detuning 
the microwaves during the interpulse period also mit- 
igates the effect of microwave leakage due to the ~25 
dB isolation of our microwave switch. Free precession is 
punctuated by the application of on-resonance irx pulses 
comprising the sequence under study. At the end of an 
experimental cycle the qubit state is rotated such that 
in the absence of dephasing, the state evolves to | j) and 
is subsequently read-out. By precisely controlling the 
length of time during which the microwaves are detuned 
during the first free-precession period (i.e., between the 
(§)x and the first ttx), we initialize the qubit in any 
state on the equatorial plane of the Bloch sphere. Con- 
trol over the final precession period is set to project the 
qubits appropriately to the dark state (in the absence of 
dephasing) . 

Measurement and operation fidelity 



given cutoff are provided a uniformly distributed random 
phase (between and 2tt) and weighted by an envelope 
function dictated by the desired noise spectrum (e.g., 
1 /uj). These components are then Fourier transformed to 
produce a time trace whose ensemble-averaged two-time 
correlation function reproduces the spectrum of interest. 
An arbitrary waveform generator outputs this time trace 
to the external frequency modulation port of a synthe- 
sizer whose beat note with a computer-controlled DDS 
produces the reference frequency for our Gunn Diode 
phase-lock loop. A peak-peak modulation V N 1 ' ax = 1 
V corresponds to a maximum relative frequency error 
of ^9xl0 -8 in a single time trace — about twice the 
linewidth of our qubit transition — and a time-averaged 
deviation of ^300 Hz at the spectral peak. We directly 
measure the noise power spectrum of our frequency- 
modulated carrier from 0.01-10 6 Hz using a phase-noise 
detection system, producing the spectra displayed in the 
insets of Fig. 3. 



Using a low-dark-count phototube we achieve a state 
measurement fidelity of ~99.5 %, limited predominantly 
by stray- light scatter from trap electrodes (Fig. 2c). We 
have microwave amplitude stability >99.9 %, and our 
pulse timing resolution of 50 ns corresponds to a system- 
atic pulse infidelity of <4xlCP 4 . A n rotation is sen- 
sitive to these inaccuracies only to second order. Care- 
ful monitoring of the resonance frequency allows us to 
maintain long-term microwave frequency stability rela- 
tive to the qubit transition of ~ 1 x 10 -9 over any mea- 
surement run, despite slow drifts and occasional sudden 
shifts in the magnetic field (hence f2). From the decay 
of Ramsey fringes we find that the integrated ambient 
noise spectrum yields a shot-to-shot magnetic field vari- 
ation 8B/B ~ 7 x 1CP 10 for 2 ms Ramsey experiments 
separated by ^200 ms. Magnetic field gradients across 
the ion array produced comparable shifts in 5B, yielding 
a 7r-pulse infidelity of <1 %. 



Nelder Mead optimization 

Experimental Nelder-Mead optimization produces 
LODD pulse sequences tailored to the user-input require- 
ments of pulse number and total free-precession period. 
Optimization begins either with PDD or the sequence 
performing best in a given noise environment. The algo- 
rithm proceeds until changes in pulse spacing converges 
to values less than the timing resolution of our pulse 
generator, generally in about 50 iterations. Optimiza- 
tion for n — 4, beginning from evenly spaced pulses in 
the ambient noise environment, yields a sequence similar 
to CPMG. We are unable to improve significantly be- 
yond CPMG, consistent with numerical simulations us- 
ing noise with a soft high-w cutoff. Optimization for a 
1 jto noise spectrum yielded a roughly 4 x improvement 
in sequence performance over CPMG at the optimization 
point. 



Noise injection 

We simulate the noise environment of other quantum 
systems by directly injecting frequency noise in our mi- 
crowave phase-lock loop. Noisy modulation of the mi- 
crowave frequency is equivalent to external field fluctua- 
tions that modulate the qubit splitting, ft, relative to a 
fixed microwave drive. By applying the noise only during 
the interpulse precession period we maintain high-fidelity 
^-rotations that otherwise would require shorter ttx- 
pulses. This technique simulates systems that apply very 
strong and narrow ixx -pulses, similar to the original pre- 
scriptions for dynamical decoupling. 

This process begins with the numerical generation of a 
desired noise spectrum; all frequency components up to a 
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